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Abstract 

Computer simulations are used to test whether a recently introduced generalization of Rosen- 
feld's excess-entropy scaling method for estimating transport coefficients in systems obeying molec- 
ular dynamics can be extended to predict long-time diffusivities in fluids of particles undergoing 
Brownian dynamics in the absence of interparticle hydrodynamic forces. Model fluids with inverse- 
power-law, Gaussian-core, and Hertzian pair interactions are considered. Within the generalized 
Rosenfeld scaling method, long-time diffusivities of ultrasoft Gaussian-core and Hertzian parti- 
cle fluids, which display anomalous trends with increasing density, are predicted (to within 20%) 
based on knowledge of interparticle interactions, excess entropy, and scaling behavior of simpler 
inverse-power-law fluids. 
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Dense suspensions of hard-sphere colloids have long served as useful experimental models 
for exploring how constraints of particle packing affect the properties of condensed phases.-"- 
Nonetheless, many technologically relevant suspensions contain high concentrations of "par- 
ticles" (e.g., micelles, microgel colloids, dendrimers, or star polymers) which interact with 
one another through considerably softer effective potentials.-""- Such systems are of funda- 
mental interest because they exhibit novel properties that cannot be readily understood 
in terms of concepts derived from idealized hard-particle systems like "excluded volume", 
"crowding", "collision rate", etc. Examples highlighted in theoretical and computer sim- 
ulation studies include reentrant and cluster crystalline phases,^"^^ as well as anomalous 
dependencies of fluid-phase thermodynamic,—"— structural,— i^^i"— and dynamioi^ii^i^SriS^i^^ 
properties on particle concentration. 

To date, much of the theoretical work on these systems has been devoted to developing and 
testing statistical mechanical approaches for predicting their structural and thermodynamic 
behavior .-i22i"— However, some of the focus is beginning to shift toward predicting dynamics. 
For example, it has recently been shown that generalized Langevin theoriea^'^'^ can quali- 
tatively (but not yet quantitatively) predict the anomalous density dependence of the long- 
time self-diffusivity observed in Brownian dynamics simulations of the Gaussian-core model. 
Likewise, mode-coupling theory has been able to capture some of the non-monotonic dy- 
namic trends displayed by fluids of particles that interact via star-polymer-like,— harmonic,— 
Hertzian,— or Gaussian-core^^i^ pair potentials. Kinetic theory has also been used to gain 
insights into the nontrivial temperature- and density-dependent trends in the long-time 
molecular dynamics of particles with bounded, penetrable-sphere interactions.— 

Another promising, albeit more heuristic, approach for predicting the dynamic proper- 
ties of soft-particle fluids is a recently proposed generalization^- of the excess-entropy scaling 
method of Rosenfeld.— >^ Here, excess entropy s*^^ refers to the difference between the en- 
tropy per particle of the fluid and that of an ideal gas of particles with the same number 
density p. Excess entropy is a negative quantity, and its magnitude characterizes the extent 
to which static interparticle correlations - present due to interparticle interactions - reduce 
the number of microstates accessible to the fluid. The qualitative expectation, corroborated 
by data from both computer simulations and experiments (see, e.g., refs. I38l-l45l). is that 
changes to a fluid that strengthen its interparticle correlations also result in slower dynamic 
relaxation processes. 
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In the generalized Rosenfeld (GR) scaling approach,— information about interparticle 
interactions is used to recast the species' long-time tracer diffusivities as dimensionless (GR- 
reduced) combinations that are, by construction, single-valued functions of s'^^ in the low p 
limit. For fluids of particles that interact solely via an inverse-power-law pair (IPL) repulsion, 
one can easily show that GR-reduced diffusivities remain single-valued functions of s'^^ for 
all values of p and temperature T. Interestingly, Krekelberg et al. have further discovered 
that GR-reduced tracer diffusivities obtained from molecular dynamics simulations of more 
complex model systems - equilibrium fluid mixtures of (additive or non-additive) hard-sphere 
or Gaussian-core particles at different values of p, T, and composition - approximately follow 
a similar excess-entropy based scaling relation.— From a practical viewpoint, this means that 
knowledge of the interparticle interactions and the thermodynamic excess entropy for these 
model fluids allows one to predict key aspects of their long-time dynamic behavior. 

Can the GR scaling method be extended to treat systems with different types of dynamics 
(e.g., with equations of motion that incorporate dynamic effects of "solvent" surrounding 
the particles)? In this work, we take a step toward addressing this question. Specifically, 
we present simulation data that tests whether the GR scaling method can be recast in a 
manner useful for predicting long-time tracer diffusivities of suspended particles undergoing 
Brownian (i.e., overdamped Langevin) dynamics in the absence of interparticle hydrody- 
namic forces. The model fluids considered here comprise particles that interact via pair 
potentials V(r) of the IPL [V(r) = e((T/r)^], Gaussian-core [V(r) = eexp{ — (r/a)^}], and 
Hertzian [V(r) = e(l — r/o")^/^] forms. Here e and a represent characteristic energy and 
length scales of these interactions, respectively, and p determines the "softness" of the IPL 
repulsion. As we discuss below, a key result from our study is that GR-reduced tracer dif- 
fusivities, sampled across a wide range of T and p for these models, approximately collapse 
onto a single curve when plotted versus s'^^. A consequence is that the tracer diffusivities 
of the ultrasoft Gaussian-core and Hertzian particle fluids, which display anomalous trends 
with increasing p,—>^ can be estimated (to within 20% relative error) - with no adjustable 
parameters- from knowledge of (1) the interparticle interactions, (2) the value of s'^^ at the 
state point of interest, and (3) the GR scaling behavior of simpler IPL fluids. 

For the case of Brownian dynamics considered here, it is convenient to express the GR- 
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reduced long-time diffusivity as=i 

(i-c/c„r^(i-c/c„) ^^-'""^^'""^i . (1) 

where B is the second virial coefficient, /3 = l/k^T, k-Q is Boltzmann's constant, and 
D2 = —{d[D/DQ]/dp)T\p=o quantifies how pair interactions at low particle concentration 
modify the tracer diffusivity D relative to the infinite dilution value Dq. As has been shown 
elsewhere,— (1 — D / Dq)'^^ = —s^^/k-Q for densities low enough that 0{p^) quantities can 
be neglected. 

Based on eq. [H it is evident that four quantities must be computed to test the GR excess- 
entropy scaling relation for state points across the T-p plane: 5, D2, D/Dq, and s'^^. The 
first two depend solely on /3 and V(r) and can be calculated from simple theoretical consider- 
ations. Specifically, B can be obtained by integrating the following expression numerically: 
B = —(1/2) f dr{exp[—f3V{r)] — 1}. Likewise, D2 can be computed via a one-dimensional 
numerical integration. 



2n prf[exp(-/3V(r))]^, , ^ 



D2= —J^ ^ Qiry dr (2) 

Eq. [2] is derived by the relaxation method, an approach outlined in detail elsewhere.—"— 
The function Q{r) is calculable from knowledge of V(r) via^ 



.d'Q dQ 

" h r — 

dr'^ dr 



2 - r.''''^^' 



2Q^r^'-^ (3) 

dr 



dr 

The physical significance of Q{r) is that it characterizes the extent to which application of a 
small external force f to a particle in the low-p fiuid would distort its pair correlation function 
with its neighbors, g{r) = exp[— /3V(r)]{l + PQ{r)r ■ f/2r}, and hence modify its effective 
mobility. One boundary condition on Q{r)— valid for all V(r) considered here - is that the 
force-induced distortion of the pair correlation function must decay at large distances from 
the central particle, Q{oo) = 0. The second boundary condition depends on the form of the 
short-range interparticle repulsion. For interactions studied previously with a hard-sphere 
repulsion at a separation r = a,—"— the following condition holds: [dQ/dr)\r=cr = —1- We 
have verified that this condition also accurately describes the behavior of particles with 
softer IPL repulsions if it is applied at an effective hard-sphere "Boltzmann"— diameter 
r = (Tb defined by V{cr^) = lOk^T. To our knowledge, this work is the first to consider D2 
for ultrasoft potentials like the Gaussian-core or Hertzian models that are bounded at zero 



interparticle separation. From eq. [3l it is apparent that the condition Q(0) = apphes for 
such models. 

To compute D/Dq, we use Brownian dynamics simulations, where particle translations 
are governed by the Langevin equation, solved in the overdamped limit. The algorithm 
we employ, wherein the position rj of each particle i G [1,A^] updates at each time step 
according to ri{t + At) = ri{t) - DqAIV Y!^_^^l3V{rij{t)) + Ar^, is detailed elsewhere.-^"^ 
The quantity Arj denotes a random displacement due to solvent collisions; its magnitude is 
taken to be Gaussian-distributed with a mean of zero and a variance of QD^At. 

The Brownian dynamics simulations track particles in a periodically-replicated cubic 
simulation cell whose volume V sets p = N/V. We use = 1000, 2000, and 3000 particles 
for the IPL, Hertzian, and Gaussian-core fluids, respectively. Interparticle potentials were 
truncated at r = lAa - 3.9a,^^ a, and 3.71a, respectively, for these systems. For simulations 
of the IPL fluids, we set Dq = Omi{a^e/mY/'^ and the At = 0.025x5, where tb = mDo/kBT 
is a characteristic Brownian time scale. For simulations of Gaussian-core and Hertzian par- 
ticle systems, we set Dq = 0.0001(cr^e/m)^/^ and At = O.lr^. Long-time tracer diffusivities 
are computed from the mean-squared displacements via the Einstein relation, D = (Ar^)/6t 
as t — 7- oo. 

We employ free-energy-based simulation techniques to evaluate s'^^. In the first step, 
we obtain the p dependence of the Helmholtz free energy at high T using grand canonical 
transition matrix Monte Carlo simulation.— In the second step, we perform a canonical tem- 
perature expanded ensemble simulation^ with a transition matrix Monte Carlo algorithnt^ 
to evaluate the change in Helmholtz free energy with T at constant p. Collectively, these 
simulations provide values for excess Helmholtz free energy and excess energy, and hence 
s^^, at a state point of interest. Interested readers can find additional details in our earlier 
papers.— "ii^ 

We now examine the simulated behavior of six IPL fluids with exponents of yU = 6, 8, 
10, 12, 18, and 36. Specifically, Figure [T^ shows how (1 — D/Dq) and its GR-reduced 
counterpart (1 — D / Dq)^^ vary with density pa^ and —s'^^/k-Q, respectively, for these fluids 
at T = e/k^. As should be expected, systems with softer interparticle repulsions (i.e., 
lower /i) display considerably faster single-particle dynamics at a given pa^. For example, 
the /i = 6 fluid must realize a density twice that of the p = 36 fluid for both systems to 
exhibit D/Dq = 0.2. In contrast, (1 — D / Dq)'^^ ioi these fluids shows only minor systematic 
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variation when plotted versus —s^^/ks, with less than 20% difference between the /z = 6 
and yU = 36 values. 

To put to the GR scaling relationship for Brownian dynamics to a more stringent test, 
we further analyze the behaviors of model fluids with ultrasoft Hertzian and Gaussian-core 
interactions (see Figures [Dd and[T]3, respectively). In particular, we explore a broad range of 
pa^ for these two systems along isotherms spanning from the moderately supercooled liquid 
(T < 0.004e/A;B) to temperatures high enough (T > COIc/Zcb) for the models to exhibit 
mean-field- like behavior Note that this data set includes state points for which D / Do 
decreases with increasing pa^ in both systems, behavior that is anomalous when compared to 
the trends exhibited by simpler atomic and colloidal fluids. The main conclusion to be drawn 
from Figures [Dd anddt is that, despite the anomalous relationship between single-particle 
dynamics and p in these systems, the GR excess entropy scaling relation looks strikingly 
similar to the quasi-universal form of the IPL fluids. In fact, the Hertzian and Gaussian-core 
data deviate by less than 20% from the least-squares flt to the (1 — D / Dq)^^ versus s'^'^ data 
from Figure [1^. 

One consequence of this quasi-universal scaling behavior is that one can estimate Brown- 
ian dynamics D/Dq data, even for nontrivial systems with ultrasoft potentials, simply from 
knowledge of their interparticle interactions, the excess entropy, and, e.g., a least-squares 
flt of the IPL data of Figure [1^. Figure |2] illustrates the quality of the predictions made by 
this approach for the Hertzian and Gaussian-core model fluids. As can be seen - with no 
adjustable parameters - the estimated values of diffusivity are accurate for low-to-moderate 
p, and remain qualitatively reliable for higher p where the anomalous trends emerge. Collec- 
tively, these results illustrate a strong apparent link between static interparticle correlations 
and long-time dynamics for model fluids with either hard or ultrasoft interactions. 

Finally, whether the approach outlined here can provide an accurate estimate to the 
"thermodynamic" component of diffusivity for systems which also have strong hydrodynamic 
interparticle contributions to dynamics is an open question which deserves attention in future 
studies. Furthermore, it will be interesting to test whether the GR scaling for diffusivity 
in Brownian dynamics also holds for systems with more complex interparticle interactions 
(e.g., with more than one characteristic length scale). 
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FIG. 1: GR-reduced long-time diffusivity, (1 — D / Dq)^^ , versus negative excess entropy —s^^/kB 
for (a) IPL, (b) Hertzian and (c) Gaussian-core fluids. The IPL fluids have exponents /i discussed 
in the text. The sohd black curve is a least-squares fit to the IPL data from Figure [1^. The dashed 
red lines indicate a difference of 20% from the IPL fit. The insets show the reduced diffusivity 
change, 1 — D/Dq; as a function of density, pa^. 
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FIG. 2: GR predictions (curves, discussed in text) and Brownian dynamics simulation data (sym- 
bols) for long-time diffusivity of the (a)Hertzian and (b) Gaussian-core fluids as a function of 
density pa^. 
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